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Abstract 

A scalar theory of gravity with a preferred reference frame is presented. It is insisted 
on the dynamics, which involves a (non-trivial) extension of Newton's second law, and 
on the new version ("v2") with isotropic space metric. We display the energy conserva- 
tion equation obtained with v2. Then the principles of the asymptotic post-Newtonian 
approximation are discussed in some detail. The results of its application to the motion 
of a small extended body in a weakly-gravitating system are given and discussed: the 
weak equivalence principle was violated in vl, due to its anisotropic space metric (as the 
standard Schwarzschild metric), but is valid with v2. 



1 Introduction 

A preferred-frame theory of gravitation, based on just a scalar field, has been previously pro- 
posed and studied by this author. It is derived from a semi- heuristic interpretation of gravity 
as due to a pressure gradient in a universal fluid or "ether" , but it consists of a definite and 
consistent set of equations; see Ref. ^ for a detailed motivation and a summary of the con- 
struction and the test of that theory. Although this theory has passed a number of tests, it has 
recently been discarded by a surprising violation of the weak equivalence principle (WEP), 
which has been found to occur for real bodies of a small but finite size |2j, whereas the WEP 
is a built-in feature of the theory for test particles. That violation can be seen, due to the use 
of a rigorous approximation method for weak gravitational fields, which the simplicity of this 
theory makes it possible to implement. 



It turns out that the theory can be saved by modifying the space metric and the equa- 
tion for the scalar field, without changing the semi-heuristic foundation nor the dynamical 
equations — which are the two characteristic features of the theory. The main aim of this con- 
tribution is to present the new version ("v2") and the asymptotic scheme of approximation 
in a self-contained way. It shall be insisted on the foundations of both the theory and the ap- 
proximation method, but many calculations shall be skipped. The reasons why the violation 
of the WEP occurs in vl and not in v2 shall be exposed. Some comparisons with general 
relativity (GR) shall also be inserted. 
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2 Presentation of the scalar ether-theory (vl &; v2) 



2.1 Preferred reference frame and space-time metric 

Space-time is assumed to be the product V = R x M, where M is the preferred reference 
body (as Newton's absolute space), endowed with an Euchdean metric g^. The equations of 
the theory are primarily written in the preferred reference frame E in which the body M is 
at rest. [A general reference frame is for us essentially a reference body N, plus a notion of 
time. As to the frame E, associated with the body M, it is endowed with the "absolute time" 
T = /c, where is the canonical projection of X £ V into R, and c is a constant — the 
limit velocity of mass points (see point 12231) ■] Thus, we shall usually formulate equations 
that are space-covariant only. The Euclidean metrics on the component spaces R and M allow 
to define a flat Lorentzian metric 7*^ : for a 4- vector U = {U^, u) (with a real number and 
u a spatial vector, i.e. an element of the tangent space TMx to M at some x G M), 

7°(U,U) = ([/0)2-30(u,u). (1) 

This metric would measure the proper time, if there were no gravity. The gravity field is a 
scalar field /3(T, x), that both has metrical effects and produces a gravity acceleration. The 
physical space-time metric 7 is related to the flat metric 7^ through the scalar field /3: 

7(U,U)=/?2([/0)2_5(u,u), (2) 

with g the physical space metric on the preferred body M, which is itself related to g^ through 
the scalar field (3 as explained in Subsect. 12.31 Equation ((JJ implies that the scalar /3 can be 
defined as 

/3^(7oo)'/^ (3) 

in any coordinates {y^) adapted to the frame E [i.e., such that each point x G M has constant 
space coordinates (y*) {i = 1,2,3)], and such that = x^ = cT. From ((21), it also follows 
that 7oi = (i = 1, 2, 3) in such coordinates. 



2.2 Dynamics 

2.2.1 The gravity acceleration in the scalar ether-theory 

An important equation of this theory is that for the gravity acceleration: 

2grad /? 

where the gradient vector with respect to the physical, Riemannian space metric g is defined 
by the space-contravariant components (gradg/3)* = with (g*-') the inverse matrix of 

matrix (gij). This equation can be obtained in two different ways. Let us outline these two 
ways. The first way is semi- heuristic: one may interpret gravity as Archimedes' thrust due 
to the gradient of the macroscopic (smoothed-out) pressure pe in an imagined fiuid called 
"micro-ether", of which the elementary particles would be local organizations (e.g. vortices). 
This gives |3] 

g = -^, (5) 

Pe 

where pe = Pe{Pe) is the macroscopic "ether density", which thus decreases towards the 
gravitational attraction. Then, in that framework, the Lorentz-Poincare interpretation of 
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special relativity (SR) |S El E] leads naturally, as summarized in Subsect. 12.31 to assume |3] 
that there are metrical effects of a gravitational field, which depend on the ratio 

a = Pe/Pf < 1, (6) 

where 

pTiT) = Sup^gM /'e(x,T) (7) 

is the ether density in remote regions that are free from gravitational field. [Precisely, this 
assumption leads to ^700 = Pe/pTi thus Eq. (jEJ is consistent with Eq. ©.] Furthermore, 
the limiting velocity c of SR is then interpreted as the "sound" velocity in the ether [H], thus 
giving 

Pe = C^Pe- (8) 

Equations © and ©, combined with ©, lead indeed to Eq. (0]). The second way to obtain 
Eq. @ is as follows: as recalled below, one may uniquely define Newton's second law in any 
relativistic theory of gravity with curved space-time, the gravitational force being m{v)g with 
m{v) = m(0)7^ the relativistic inertial mass and g a theory-dependent gravity acceleration. 
(t> = gi(v, v)^/^ is the modulus of the velocity v, and 7^ is the Lorentz factor.) Now, /? being 
defined by Eq. (jSJ, postulating Eq. @ for vector g is equivalent, under natural requirements, 
to asking geodesic motion in the case of a static gravitational field jlj . 



2.2.2 Newton's second law in a curved space-time 

To use a gravity acceleration, we need to define dynamics of a test particle by extending 
the special-relativistic form of Newton's second law to the case with gravitational force in a 
curved space-time — instead of postulating that free test particles follow space-time geodesies, 
as is done in GR. This can be done for any theory with Lorentzian space-time metric. The 
gravitational force has just been introduced, and the momentum must of course be 

P=m{v)v, (9) 

where the velocity v has to be defined in terms of the "local time" tx, measured by a clock 
at the spatial position x = (y*) that is fixed in the frame F considered, and that momentarily 
coincides with the position of the test particle — thus 

V = dx/dtx- (10) 

{In any theory, the frame F is fixed by the choice of the space-time coordinate system (y^), 
though many coordinate systems correspond to the same frame. The corresponding body N 
is described by the space coordinates (y*), which can be subjected to purely spatial changes 
yH — ^t^yj-j [TllSlin].} For the scalar ether-theory, we consider only the preferred reference 
frame, thus F = E, and we have from 

dWdT = /3(r,x). (11) 

However, a synchronized local time may also be defined for any trajectory ^ 1— > X(^) = {y^{(,)) 
in the arbitrary reference frame F, in any theory with curved space-time metric {see Landau 
& Lifshitz |B1, Eqs. (84,14) and (88,10), and Cattaneo 0}: 

dtx _ Vloo f dy° 7o» dy'\ 

c V dc 700 de y ■ ^ ^ 
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To complete the definition of Newton's second law, we have to define the time-derivative of a 
spatial vector (the momentum P of the test particle in the frame F). This is not trivial in the 
general case of a variable gravitational field, which means a time-dependence of the spatial 
metric g in the frame F. We have found [71 that a unique definition may be given for the 
time-derivative Umv/D^ of a space vector w depending on a parameter £^ along a trajectory in 
a space (manifold) N endowed with a metric g = g^ that varies with the parameter ^, under 
compelling requirements that include the validity of Leibniz' differentiation rule for a scalar 
product. In coordinates, this unique definition is as follows: 

with r*-^ the Christoffel symbols of metric g. Our extension of Newton's second law is hence: 

DP 

m{v)g = —, (14) 

thus the parameter ^ in (|13j) is then the synchronized local time, given in general by Eq. (|12j) . 
and in the investigated theory by Eq. It turns out that Einstein's motion following the 

geodesic lines of 7 can be put in the form ^14^ : if the coordinate system can be chosen such 
that 7oi = {i = 1,2,3) (which is true for a class of possible reference frames), the gravity 
acceleration has to be (Tj 

2gnKy? 1 

ggcod = -c — -t.v, (15) 

where (3 = ^/joo and g is the space metric in the reference frame fixed by the coordinate 
system. This is invariant under the coordinate transformations [7] that leave the reference 
frame unchanged. 



2.2.3 The limit velocity and the case with photons. The geodesic motion 

It follows from Eq. (|14() that a "mass point" (a test particle with positive rest-mass) cannot 
reach the velocity c. Hence, that velocity is indeed an upper limit for mass points, as in 
special relativity. For the investigated theory, this limit applies a priori in the ether frame E, 
however, because the equations for the gravity acceleration (jl} and for the metric ((2)) apply 
in that frame. For a photon, the same law is used, substituting in Eqs. © and ()14() the 
mass content E/c^ of the energy E = hv for the inertial mass m{v), the frequency v being 
evaluated with the local time. As applied to photons, our extension of Newton's second law 
implies that they keep the constant velocity c (again in the frame E) {see Eq. (29) in Ref. 
jl()|}. Now, when going to an arbitrary reference frame, we transform the space-time metric 
[defined in the frame E by Eq. 0] as a space-time tensor, of course. This actually means 
that we are using the standard (Poincare-Einstein) synchronization convention of clocks {see 
Ref. jSl, Sect. 84}. Thus the time-like, light-like, or space-like character of a world line 
element is, of course, invariant. In other words, photons will also run at the velocity c in any 
physically possible reference frame, but, in frames other than the ether frame, this has partly 
a conventional character: only the hack- and- forth velocity of light is physically meaningful 
[H El As already mentioned, it also follows from Newton's second law H14() with the 
gravity acceleration Q) that, in the case of a static gravitational field (which, in that theory, 
is characterized by condition [3^ = 0), test particles follow geodesic lines of metric 7. This is 
true for mass points [7j [as can be seen immediately by comparing Eqs. Q) and (fT3)) ] and for 
photons as well jl2j . 
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2.2.4 Continuum dynamics and Maxwell equations 

To describe the motion of a continuous medium (be it a fluid or an electromagnetic field, 
for instance), we begin by investigating the case of a dust, i.e., a continuum made of non- 
interacting mass particles. For a dust, Newton's second law (|14|) with the gravity acceleration 
Q may be applied pointwise and implies |13j that the material energy-momentum tensor T 
verifies the following equation: 

TH;. = ft/., (16) 

where 

[Indices are raised and lowered with metric 7, unless mentioned otherwise. Semicolon means 
covariant differentiation using the Christoffel connection associated with metric 7. Note that 
in GR, in contrast, we have 6^ = in ((1^ .] The universality of gravity and the mass-energy 
equivalence mean that the same equation must hold true, with the same definition (|17j) . 
for any continuous medium. In particular, Eq. H16|) may be applied to the situation with 
both a charged medium and the corresponding electromagnetic (e.m.) field. It is in that way 
that the gravitational modification of the Maxwell equations is got in the scalar ether-theory 
|14| . This modification is fully consistent [14 with photon dynamics as defined [with E/c^ 
in the place of m{v)] by Eq. (|14() with (jlj). The wave equivalent of free light-like particles 
are the null e.m. fields in vacuo. The consistency means that the trajectories of the energy 
flux of such fields coincide with the trajectories of free light-like particles Jl]. This entails, 
in particular, the prediction that, for null e.m. fields in vacuo, the energy flux has velocity 
c. When the e.m. field is not a null field, it is not immediate how to define the velocity of 
the energy flux. It is worth mentioning that Nimtz and his team claim to have measured 
genuinely superluminal signal velocities in "photonic tunneling" experiments |15j . 

2.3 Specific form for tfie spatial metric 

The foregoing dynamics holds independently of any specific form for the spatial metric g in 
the frame E, that enters Eq. (jJJ. If one follows the semi-heuristic analysis of gravity as due to 
the macroscopic part of the pressure force in the ether, outlined at the beginning of Subsect. 
12.21 then one expects a gravitational "rod contraction" in the same ratio (j = ^/too a-s the 
gravitational "clock retardation" In short, if one analyses Special Relativity within 

that concept of a barotropic ether, one finds that an observer moving at velocity u with respect 
to E should consider that, due to the Lorentz contraction, the ether density is smaller for her 
than for a fixed observer, p'^ = Pe/'Ju- Thus, the Lorentz contraction of her meters, as well 
as the Larmor dilation of the period of her clock, depend on precisely the ratio of the ether 
density in her moving frame to that in the ether frame. A gravitational field, characterized 
(according to this semi-heuristic analysis) by a decrease of the density pe in the direction of 
the gravity acceleration, should hence also lead to rod contraction and clock retardation in 
the same ratio /? = ^/^/oo = Pel ■ I-e., the physical (g) and Euclidean (g^) space metrics 
on the preferred body M should be related together by a dilation of physically-measured 
distances, with respect to Euclidean distances (due to the contraction of objects, including 
length standards). But that dilation can occur: 

• either in a locally-isotropic way, which means postulating the following relationship 
between the two spatial metrics 0|: 

g = rV (v2). (18) 
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• or, as in SR, in just one direction — which then can only be that of the gravity acceleration 
g. This was the assumption set in the first version ("vl") of the scalar ether-theory, 
which passed a number of tests but which seems to be discarded by the violation 

of the weak equivalence principle discussed in Sect. |S1 below. 

Therefore, the author is now investigating in detail the new version ("v2"), based on Eq. H18|) . 



2.4 The question of energy conservation 

The conservation of energy is at the very top of the hierarchy of the physical laws. Unfortu- 
nately, there is no exact local conservation law for energy in GR. This is due to the fact that 
energy is not a generally-covariant concept. Essentially, the equation for continuum dynamics 
in GR, r^.j^ = 0, may be rewritten as the cancellation of the 4-divergence of a tensor-like 
object with respect to a flat metric — i.e., it may be rewritten as an exact local conservation 
law. But the equation obtained thus is covariant only under linear coordinate changes, and, 
moreover, the corresponding "gravitational energy-momentum pseudo-tensor" is not unique 
(see e.g. Stephani jEl; see Ref. ^H] for this author's opinion and some other references). 

Since the theory investigated in this paper is unshamedly a preferred-frame theory, one 
should be able to ask for a local energy conservation in this theory. To investigate this 
question, we rewrite the time component of the equation for continuum dynamics (|16j) in 
terms of usual derivatives: 

(^Tj) . + iV^T^)^, = V^^oT'^, (19) 

where 

7 = det(7^^) = -P^det{gij). (20) 

On the l.h.s. of Eq. (|19|) . we have the energy balance of matter, whereas on the r.h.s. we 
have a source term depending on the matter energy density (in the preferred frame), T^^, 
and on the scalar gravitational field, (3 = ^/joo- Note that, in the sought-for equation for 
the scalar gravitational field, /? might be replaced by some function of it, for convenience. 
This function and the equation for the scalar field are not fully constrained by the heuristic 
principles of the ether theory Another constraint will thus be that, by virtue of this scalar 
field equation, one should be able to transform the r.h.s. of H19|) to a 4-divergence, in order to 
get a conservation equation for a (material plus gravitational) energy. Due to Eq. H2U() . this 
depends on the form assumed for g. 



2.5 The equation for the scalar field and the energy conservation 

The heuristic interpretation of gravity as Archimedes' thrust in a compressible fluid leads |S| 
to postulate the following general form for the equation governing the scalar gravitational 

Pe + (time derivatives) = — —apeF{(3), F{(3) ^ 1 as /? — > 1, (21) 

where Ag is the Laplace-Beltrami operator, i.e., the Laplacian associated with the Rieman- 
nian metric g. The "time derivatives" term should be such that, in the appropriate limit, 
involving the condition /? — > 1 uniformly, the operator on the l.h.s. should become equivalent 
to the usual (d'Alembert) wave operator. [As is seen on Eqs. Q and (fTH)) . (3^1 means that 
the space-time metric 7 tends towards the flat one "f^.] On the r.h.s. of (|21|) . a denotes the 
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relevant mass-energy density, which has to be defined in terms of the energy-momentum ten- 
sor T. Obviously, this nonlinear wave equation is not fixed. In the first version of the theory, 
a definite equation had been postulated involving the particular condition F{(3) = 1. 
Moreover, it had been found that the assumed equation could be rewritten in terms of 
the flat metric 7^ and the field f = 0^, and did imply a local energy conservation. 

With the space metric (|18|) assumed in v2, one finds that the term Ag pe in Eq. (|21|) can 
be reexpressed in terms of the Euclidean metric, namely as — AgOip, with = —Log/?. 
And one finds that the r.h.s. of Eq. ^ can be rewritten as —xj^flT^^. Moreover, if one 
assumes the following linear wave equation for the scalar field [corresponding to F{f3) = P'^ in 

Eq. mv- 

ni; = ^^,0,0 -AiP = (22) 

with 

a = r°°, (23) 
then one obtains a local conservation equation for the energy: 

doiSra + Sg) + div(*m + *g) = 0. (24) 

In Eq. (|24() . the material and gravitational energy densities are given (in mass units) by: 

,2 



TOO, e - 



8ttG 



^2 +(gradV')' , (25) 



and the corresponding fluxes are: 



^ (TO^), *g ^ (V,o grad^) . (26) 

(Henceforth, the operators grad, div and A shall be relative to the Euclidean metric g^.) We do 
assume Eq. (|22j) as the equation for the scalar gravitational field in v2. This equation and the 
local energy conservation H24|) are invariant under arbitrary changes of spatial coordinates, but 
need that the time coordinate is the preferred one, = cT. Thus, these are preferred-frame 
equations. The energy conservation ()24() substitutes for the mass conservation: independently 
of the assumed form for the space metric, and thus in v2 as well as in vl, one may indeed 
show that the dynamical equations (jlGj) imply a reversible creation / destruction of matter in 
a variable gravitational field [EI- However, the theory says that mass is extremely close to 
be conserved in usual situations fH^ . 



3 Asymptotic post-Newtonian approximation (PNA) 

3.1 Principle of an asymptotic scheme for the field equations of gravitation 

Very accurate experimental comparisons are aimed at in the theory of gravitation, hence 
one should use a very clean approximation scheme. The purpose of the PNA is to obtain 
asymptotic expansions of the fields as functions of a small parameter A (the field-strength) . A 
clean way to obtain this, in accordance with the general principles of asymptotic analysis, is 
to deduce a family of systems (fields), (S'*'), from the data of the system of interest, S, and this 
in a general-enough case. A system is formally defined as the solution of a boundary- value 
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problem for the set of field equations. Since gravitation propagates with a finite velocity in 
relativistic theories of gravitation, the relevant boundary-value problem is the initial-value 
problem. Hence one has to deduce a family of initial data from the initial data valid for 
the system of interest, S (e.g. the solar system). Then one may expect that the solution 
fields do admit asymptotic expansions — usually Taylor expansions — as A — > 0. By inserting 
these expansions into the field equations, one gets each of the latter ones "split" to successive 
equations corresponding to the increasing orders in A in the Taylor expansions. The expanded 
equations, obtained thus, are much simpler than the original equations. The system of interest, 
S, is assumed to correspond to a small value Aq <C 1 of A, so that one may use the expansions 
and the expanded equations for S itself — which, of course, is the real aim of the expansion 
method. 

3.2 Physical framework for the Newtonian hmit 

Naturally, the family (S'*') of systems should be relevant to the physical situation applying to 
the system S, which situation justifies the use of asymptotic analysis. In the case of the PNA, 
the situation is that the gravitational field is weak, and accordingly Newtonian gravity (NG) 
is supposed to give a very accurate first approximation to the exact behaviour. In less vague 
terms: 

• i) the physical space metric is close to being Euclidean and the local physical time is 
close to flowing uniformly (equivalently: the physical space-time metric 7 is close to the 
flat metric 7*^); 

• ii) the matter fields are nearly equal to solution fields of the Newtonian equations; 

• iii) the Newtonian potential [/n (or rather some equivalent V of U-^, defined internally 
to the relativistic theory considered) is small. 

To make conditions ii) and iii) precise, it is necessary, as a preliminary step, to define a 
purely Newtonian situation, in which C/n itself would be small. In other words, one is led 
to ask whether a weak-field limit of NG itself could be defined, by considering a family (S^) 
of Newtonian systems. One then finds |2n[ I2T] that there is a similarity transformation in 
NG, which is exactly appropriate to that very purpose. If p^^^ (pressure), u^^) (velocity), 
p^^^ = F^^\p^^^) (density), and U^'' obey the Euler-Newton equations (the field equations for 
a perfect fiuid in NG), then, for any A > 0, the fields 

pW(x,r) = AV'^x, VA T), pW(x,T) = Ap(^)(x, VA T), (27) 

4^)(x,r) = AC/jf)(x,VAr), uW(x,T) = ^u«(x,VAr), (28) 

also obey these equations — provided the state equation for system is F^^^ (p^"^) ) = AF^-*^) ( A^^p^"^) ) . 
As A — > 0, the potential and the density in the bodies decrease like A (while the bodies keep 
the same size), the velocities decrease like \/A, and accordingly the time scale increases like 
1/\/A. Thus, this similarity transformation defines indeed the weak-field limit in NG itself. 
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3.3 Definition of the family of gravitating systems 

Therefore, the idea is to define the initial data of system S^, in the relativistic theory consid- 
ered, by applying the similarity transformation H27j) - ()28|) to the initial data for system S. ^ If 
conditions i) to iii) above are asymptotically satisfied as A — > 0, this will undoubtedly define 
a relevant family (S^) of systems. However, one has to find an equivalent of the Newtonian 
potential in the theory considered. In the scalar ether-theory, condition i) means exactly 
that the scalar field ^ = —Log ^/too must tend towards (for v2. For vl it meant that the 
field / = 7oo had to tend towards 1 In the asymptotic framework, and in view of the 

Newtonian similarity transformation, condition ii) needs the validity of the following one: 

• ii') The matter fields in system have the same order in A as in the transformation 

(Ulli-dll. 

If that condition applies to a family (S^), then, by changing the mass and time units for 
system S'^ in this way: [M];s^ = A[M] and [T];^ = [T]/\/A, the matter fields p^'^\ p^^\ and 
u^'^\ become ord(A''). Moreover, in these units, we have A oc With the help of these 

units, one finds easily that, if the field equation ()22|) is to be satisfied by each member of the 
family, then the field ip^"^^ must be order 1 in A, like the Newtonian potential C/n. In this 
sense, ip is an equivalent of U-^ in v2. Since ip (unlike U-^) is dimensionless, the most natural 
field-strength parameter is simply 

A = SupxeM V'(x) (29) 
(at the initial time T = 0, say). The field 

V = c^V (30) 

will also be ord(A) (in fixed units) if ip is ord(A), in addition it satisfies the wave equation 
with the same r.h.s. (in the Newtonian limit where u ~ p) as Poisson's equation of NG. Since 
the retardation effects should become negligible in the Newtonian limit, V^"^^ will be really 
equivalent to as A — > 0. Finally, u = ord(\/A) (as demanded by condition ii') means 
a A~^/^ dependence of the characteristic time; if one combines this with the fact that V is 
ord(A), one expects the asymptotic validity of Eq. (|28|) i for V : 

y(^)(x,T) ~ Ay(^)(x, VAT) as A ^ 0. (31) 

In uniform conditions, one may differentiate (|3H) with respect to T: 

dTV^^\x,T) ~ X^^^dTV^^\x,VXT) as A ^ 0. (32) 

Therefore, in order to satisfy conditions i), ii') and iii), it is most natural to define the family 
(S^) by the following family of initial conditions: 

pW(x) = A2p«(x), p*W(x)=Ap*W(x), (33) 

uW(x) =^/Au«(x), (34) 

^ Or more precisely, since S is assumed to correspond to a small value Ao ^ 1, by applying first the 
transformation to go from Ao to A = 1, and second from A = 1 to the arbitrary value A. This amounts to 
substituting ^ = A/Ao for A, and p^^°\ etc., for p^^\ etc. |21|. 



9 



yW(x) = Ay(i)(x), dTV'^^\x) = X^/^dTV^^\x). (35) 

{p* is the proper rest-mass density.) We may confidently expect that, in a significant interval 
of the "dynamical time" |2Uj r = T\/A, the solution fields will remain of the same order in A as 
in this initial data, hence conditions i), ii') and iii) will indeed be satisfied in the asymptotic 
sense, in that interval. 

3.4 Implementation of the method in the scalar theory and in GR 

Using the varying units [M]x = A[M] and [T];^ = [T]/\/A, the matter fields and the scalar 
gravitational field V become ord(A^), while the small parameter A oc This makes the 

expansions straightforward. All fields (of course) depend on A and are Taylor-expanded {at 
least up to some order n^aaxi beyond which other expansion functions might have to be used, 
as done in GR (22, ,22, } : 

n 

p = ^p^A'= + 0(A"+i), etc. (andA = l/c2) (36) 
A;=0 

Inserting these expansions into the exact equations, each of the latter ones splits to n + 1 
exact equations (it is just coefficient identification and remains true if the expansion basis 
is larger). Hence, the numbers of the independent equations and the independent unknowns 
remain equal. The zero-order equations (corresponding to /c = 0) turn out to be the Euler- 
Newton equations of Newtonian gravity. Thus, the Newtonian limit is satisfied, in particular 
it is so for condition ii). The {k = l)-equations give the first PN (IPN) correction, which is 
linear in the IPN fields, like pi. The details can be found in Ref. where the method is 
applied to vl. The modification for v2 is straightforward. 

In GR, a method similar to that one which we have just presented, has been introduced in 
a particular case by Futamase & Schutz |^ (although the present author came to this method 
independently [2B). Namely, their initial condition for the matter fields was the same as here, 
Eqs. (|33() - H34|) . but their initial condition for the space metric and its time-derivative, playing 
the role of Eq. (|35jl here, was their Eqs. (3. 13)4-5, which are in our notation: 

V^f^{^) = 6^J, (V^7^J|q(x) =0 (l<i<3, 1<J<3). (37) 

This condition is verified when the spatial metric has the "conformally-Euclidean" (or isotropic) 
form 

gij = V-7 ^ij, (38) 

in fact it is just the condition that one would impose if one would wish to have the space metric 
(|38|) . without imposing any a priori restriction on the factor ^—7. This initial condition does 
not directly impose that the gravitational field is weak, since H37() . in contrast with ()28() i and 
does not depend on the small parameter A. This seems a priori surprising. Thus, in the 
harmonic gauge, assuming the metric H38|) together with the Newtonian-limit conditions for 
the matter fields, (|33|) and (|34j) . constrains the solution of the Einstein equations to be a weak 
field. In the scalar theory also, we used the field equation to get condition but things 
look less straightforward in GR. In the framework of GR, it is legitimate to assume the initial 
condition (|S7j) . but it looks quite particular: there are many, very different, possibilities, and 
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one would like to know what would happen if he would choose very different initial conditions 
for the space metric, e.g. ones corresponding to an "anisotropic" space metric. (We shall 
come back to this question in Subsect. 15.41 ) ^ Finally, Futamase & Schutz [201 derived the 
local equations in a form which is not explicit enough for further use, and derived no "global" 
equations i.e. no equations of motion for the mass centers of extended bodies. 

3.5 Comparison with the standard PNA scheme used in GR 

An extremely good agreement is found between astronomical observations and calculations 
based on GR ;25l. However, such calculations use the standard PNA scheme, which remains 
essentially identical to that one which was proposed by Fock j26j and Chandrasekhar |27j. In 
the standard scheme, no family of systems is envisaged; is formally considered as a small 
parameter; and the matter fields p, p, u are not expanded. Since, in an asymptotic framework, 
the matter fields do depend on the small parameter, it follows that the splitting of each exact 
equation into n + 1 equations would not be justified in the standard scheme. The problem 
is that this splitting is needed at least for the equations for the gravitational field ^1 1^ . 
Detailed arguments, and references to more recent literature, can be found there |17| I21j. 

4 Equations of motion of the mass centers 

The mass centers of the celestial bodies are defined [2H] as local barycenters of the rest-mass 
density. Therefore, to get the PN equations of motion of the mass centers (EMMC's), the 
local PN equations of motion [which are just the expansion of Eq. H16() . specialized to a perfect 
fluid] are integrated in the domain Dq occupied by body (a) . One finds that the internal 
structure of the bodies (e.g. the density and velocity profiles inside them) influences the motion 
already from the first PNA [20] • This follows naturally 29 from using the "asymptotic" 
method of PNA and, in the author's opinion, should hold true for GR if an asymptotic PNA 
was used there. 

5 Point-particle limit of the EMMC's and the WEP 

Once it had been recognized [2^1 that the internal structure of the gravitating bodies does 
influence their motion, it was important to check whether or not this influence cancels in the 
limit of a very small body. If it does not, the weak equivalence principle (WEP) is violated. 

^ When trying to implement an asymptotic PNA scheme for GR, there appear two basic difiiculties as 
compared with the present theory, which is much simpler than GR. One arises from the fact that the grav- 
itational field is then a space-time tensor, thus has much more freedom than in Newton's theory, so that a 
priori the Newtonian limit might be defined in different ways. The other difficulty is due to the fact that the 
(/lO) components of the Einstein equations C^" — kT'^" are not evolution equations, and instead impose four 
nonlinear constraints on the initial data for the metric |24|. Depending on the gauge condition adopted, it 
may be more or less difficult to account for these "constraint equations" in the definition of the initial data. In 
Ref. |2U| . the harmonic gauge was used. In that well-known gauge, the Einstein equations, thus including the 
constraint equations, take a relatively simple form, but still they remain nonlinear. As a matter of fact, the 
initial conditions for the -^—7 7'''' components (which conditions have to hold as a consequence of and 
the Einstein equations) were not solved explicitly in Ref. |20| . 
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5.1 Framework for the point-particle limit 

We consider a system of N bodies (1), (A^), that move under the IPN gravitational field 
produced by them all. (Since, as explained in Subsect. 13.41 above, the asymptotic IPN 
approximation leads to split each of the exact local equations to two equations, the local 
equations of the first PNA make a closed exact system j^. Hence, assuming that the physical 
system under study is described accurately enough by the first PNA, we do not need to 
consider the weak-field parameter A.) The size of one of the bodies, say (1), is a small 
parameter ^, with ^ — > 0, thus a family (S^) of IPN systems is considered The initial data 
is independent of ^ apart from the size of the small body. One just has to find the main terms 
in the expansions as ^ — > of the integrals PB_ that enter the general form of the EMMC's. 
This is a bit long and technical but not really difficult 

5.2 Result of the point-particle limit 

The point-particle limit of the (IPN) EMMC's of vl was taken 0. It was found that a 
structure-dependent part of the acceleration. As, remains at the point-particle limit ^ — > 0. 
Moreover, in the static spherically symmetric (SSS) case, the limit value of the acceleration 
as ^ — > differs from the acceleration of a test particle just by A5: 

lim A^ = As + (IPN acceleration of a test particle in the SSS field). (39) 

Thus, the Weak Equivalence Principle is violated in vl. This was a bad surprise, because the 
validity of the WEP for test particles and even for a "test continuum" was built in the very 
construction of the theory, see Eqs. (|14|) and (|16jl above. Moreover, this violation seemed very 
likely to kill vl, since, for instance, | A5I is of the order of 10^ times the anomalous residual 
acceleration found |Hn[l^ for the Pioneer 10 spacecraft... 

Therefore, the reason for the violation was carefully analysed |2j (see Subsect. 15.31 below 
for a summary), and v2 was constructed in order to avoid this violation. And indeed the 
violation is found to not take place (at the IPN approximation) in v2: we have A5 = in the 
general case for v2, this has been checked. Moreover, Eq. (|39|) still holds (thus with A^ = 0) 
for v2. I.e., for v2, the acceleration of an extended body (1) in the SSS field produced by 
a massive body (2) coincides, as closely as desired, with the acceleration of a test particle in 
the same SSS field, provided that body (1) is small enough. 

5.3 Reasons for the WEP violation in vl and for no violation in v2 

The general reason which makes it a priori possible that the WEP might be violated for ex- 
tended bodies is simply the non-linearity of the theory! Indeed, any kind of material object 
must contribute to the gravitational field. In Newton's theory, the gravitational attractions 
between the particles of a body sum to zero (and the same is true for the corresponding 
torques), due to actio-reactio equality. (It can also be checked on the continuous form of 
Newton's attraction law, based on the Poisson equation, that the resultant of the "self field" 
created by a body is zero.) However, the actio-reactio equality cannot even be formulated in 
a nonlinear theory, because one cannot define what is the "gravitational attraction exerted by 
an object over another object". Then, the body contributes nonlinearly to the gravitational 
field and it cannot be guaranteed in advance that the body's contribution to its own acceler- 
ation will decay with the body's size. 
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Moreover, the influence of the structure of the bodies on the gravitational field, hence 
on the motion, is to be generally expected in any relativistic theory, due to the mass-energy 
equivalence. The asymptotic PNA separates the equations of motion of the different orders. 
This makes the structure influence occur explicitly in the general equations for the IPN cor- 
rections to the motion of the mass centers, because these equations involve several integrals of 
the Newtonian matter fields [2H1- Again, it is not a priori obvious that this structure infiuence 
must cancel in the point-particle limit. 

But the specific reason that makes the WEP violation actually occur in vl is the presence 
in vl's PN spatial metric g^^i-^ of the spatial derivatives U^i of the Newtonian potential. Indeed, 
the local equation of motion depends of course on the Christoffel symbols of the metric, thus, 
in that case, it depends on the second derivatives U^ij. The equation of motion of the mass 
center a of body (1), being got by integration of the local equation inside the volume of (1), 
depends on an integral involving the C^,jj's. Now the self part of U^ij [i.e., the contribution 
of (1) to U^ij] does not depend on the size ^ of (1), hence does not evanesce with ^. This is 
the reason why the acceleration a depends on the structure of the small body (1) itself, even 
in the hmit C ^ 0. 

Since, in v2, the spatial metric is isotropic [Eq. 1)181) ]. its PN expansion involves the 
Newtonian potential U (precisely: U = Vq, which obeys indeed the Poisson equation), but it 
does not involve its spatial derivatives U^i. Therefore, the WEP should not be violated any 
more in the point-particle limit of the IPN EMMC's in v2. However, the scalar field equation 
and the metric being different from vl, this had to be checked carefully. This has been done, 
and indeed the WEP applies (at the IPN approximation) to v2. 

5.4 The case with general relativity 

Now the C/,i's, which are responsible for the WEP violation in vl of the scalar ether-theory, 
are also there in the PN approximation to the standard form of Schwarzschild's metric. The 
latter is also spatially anisotropic, in just the same way as is the metric of vl (i.e., the direc- 
tion g, here the radial direction, is preferred). In fact, the standard Schwarzschild metric is 
the unique solution of vl for the SSS case (under the usual condition that the behaviour is 
Newtonian at infinity) — see e.g. Ref. ^0], Subsect. 2.1. 

Hence, it is natural to ask the question: what about GR? In the usually-adopted harmonic 
gauge, the standard IPN approximation of the spatial metric is isotropic and contains U but 
not the U/s, hence no WEP violation does occur. Now in the frame of an asymptotic PNA 
scheme, imposing an initial condition corresponding to an isotropic metric is consistent with 
the harmonic gauge [201) but it should be tried to impose different conditions (Subsect. Il-{.4j) . 
Moreover, gauges do exist, in which the standard Schwarzschild solution is the unique solution 
of the SSS case with Newtonian behaviour at infinity (e.g. one is shown in Ref. Note 15; 
a more convenient one is also known to the author). In a such "Schwarzschild- like gauge", 
the spatial metric of the general case would be expected to be "anisotropic", and its PN 
approximation would be expected to involve derivatives U^i. Hence, it seems quite likely that 
the WEP could also be violated in GR with a such gauge. In GR, then, the validity of 
the WEP for a small extended body at the IPN approximation might depend on the gauge 
condition. 
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6 Conclusion 



A scalar theory with a preferred reference frame has been summarized. A new version ("v2") 
has been presented. In vl, the space metric was anisotropic, whereas in v2 it is isotropic. 
The complete metric of v2, Eq. with H18|). is an often- met one [5^2. 33, 34 . The formal 
specificity of the present theory (independently of its motivation and its heuristic) is in its 
revendicated preferred-frame character (Subsect. 12. its dynamics (Subsect. 12.2(1 . and its 
scalar field equation implying a local energy conservation (Subsect. 12.5(1 . 

To test that theory in celestial mechanics, an "asymptotic" PN scheme, as in applied 
mathematics, was developed. The resulting equations of motion for a self-gravitating system 
of extended bodies include internal- structure effects. For vl only, the internal-structure influ- 
ence subsists at the point-particle limit. This is a severe violation of the WEP, although the 
WEP is satisfied for test particles. v2 was constructed precisely to avoid this violation — and 
indeed it does avoid it. 

Using an "asymptotic" approximation scheme, the same WEP violation as in vl might 
occur in GR also, e.g. in a gauge where the SSS solution is the standard Schwarzschild metric, 
and/or with an initial condition corresponding to an anisotropic space metric — but it would 
be more difficult to check it than in the present theory, due to the greater complexity of GR. 
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